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Transverse Momentum Distribution and Nuclear Modification Factor of 
Charged Particles in p-Pb Collisions at y^s^ = 5.02 TeV 



The transverse momentum (pi) distribution of primary charged particles is measured in non-single- 
diffractive p-Pb collisions at ^/Imn — 5.02 TeV with the ALICE detector at the LHC. The pj spectra 
measured near central rapidity in the range 0.5 < pj < 20 GeV/c exhibit a weak pseudorapidity de- 
pendence. The nuclear modification factor Rpp\, is consistent with unity for pj above 2 GeV/c. This 
measurement indicates that the strong suppression of hadron production at high pj observed in Pb-Pb 
collisions at the LHC is not due to an initial-state effect. The measurement is compared to theoretical 
calculations. 
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Measurements of particle production in proton-nucleus collisions at high energies allow the study of 
fundamental properties of Quantum Chromodynamics (QCD) at low parton fractional momentum x and 
high gluon densities (see 1 1 ] for a recent review). They also provide a reference measurement for the 
studies of deconfined matter created in nucleus-nucleus collisions li2l . 

Parton energy loss in hot QCD matter is expected to lead to a modification of energetic jets in this medium 
(jet quenching) [T]. Originating from energetic partons produced in initial hard collisions, hadrons at high 
transverse momentum pi are an important observable for the study of deconfined matter. Experiments 
at RHIC have shown |4 , 5 1 that the production of charged hadrons at high pi in Au-Au collisions is sup- 
pressed compared to the expectation from an independent superposition of nucleon-nucleon collisions 
(binary collision scaling). 

By colliding Pb nuclei at the LHC it was shown that the production of charged hadrons in central 
collisions at a center-of-mass (cms) collision energy per nucleon pair ^snn = 2.76 TeV shows a stronger 
suppression than at RHIC, indicating a state of QCD matter with an even higher energy density. At the 
LHC, the suppression remains substantial up to 100 GeV/c QH and is also seen in reconstructed jets Q. 
A p-Pb control experiment is needed to establish whether the initial state of the colliding nuclei plays a 
role in the observed suppression of hadron production at high-/7T in Pb-Pb collisions. In addition, p-Pb 
data should also provide tests of models that describe QCD matter at high gluon density, giving insight 
into phenomena such as parton shadowing or gluon saturation HI. 

In this letter, we present a measurement of the pi distributions of charged particles in p-Pb collisions 
at -y/^NN = 5.02 TeV. The data were recorded with the ALICE detector |10| during a short LHC p-Pb 
run performed in September 2012 in preparation for the main run scheduled at the beginning of 2013. 
Each beam contained 13 bunches; 8 pairs of bunches were colliding in the ALICE interaction region, 
providing a luminosity of about 8 x 10^^ cm~^s~^ The interaction region had an r.m.s. width of 6.3 cm 
in the longitudinal direction and of about 60 /xm in the transverse directions. 

The trigger required a signal in either of two arrays of 32 scintillator tiles each, covering full azimuth and 
2.8 < Tjiab < 5.1 (VZERO-A) and -3.7 < Tjiab < -1.7 (VZERO-C), respectively. The pseudorapidity in 
the detector reference frame, Tjiab = — ln[tan(0/2)], with Q the polar angle between the charged particle 
and the beam axis, is defined such that the proton beam has negative Tjiab. This configuration led to 
a trigger rate of about 200 Hz, with a hadronic collision rate of about 150 Hz. The efficiency of the 
VZERO trigger was estimated from a control sample of events triggered by signals from two Zero Degree 
Calorimeters (ZDC) positioned symmetrically at 112.5 m from the interaction point, with an energy 
resolution of about 20% for single neutrons of a few TeV. 

The offline event selection is identical to that used for the analysis of charged-particle pseudorapidity 
density (dA'^^h/dT]iab) reported in 111]. A signal is required in both VZERO-A and VZERO-C. Beam- 
gas and other machine-induced background events with deposited energy above the thresholds in the 
VZERO or ZDC detectors are suppressed by requiring the signal timing to be compatible with that 
of a nominal p-Pb interaction. The remaining background after these requirements is estimated from 
triggers on non-colliding bunches, and found to be negligible. The resulting sample of events consists 
of non-single-diffractive (NSD) collisions as well as single-diffractive and electromagnetic interactions. 
The efficiency of the trigger and offline event selection for the different interactions is estimated using 
a combination of event generators, see [11] for details. An efficiency of 99.2% for NSD collisions is 
estimated, with a negligible contamination from single-diffractive and electromagnetic interactions. The 
number of events used for the analysis is 1.7 x 10^. 

The primary vertex position is determined with tracks reconstructed in the Inner Tracking System (ITS) 
and the Time Projection Chamber (TPC) by using a 'x}' minimization procedure described in [8|. The 
event vertex reconstruction algorithm is fully efficient for events with at least one track in the acceptance, 
|T7iab| < 1.4 (when the center of the interaction region is included as an additional constraint). An event 
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is accepted if the coordinate of the reconstructed vertex measured along the beam direction is within 
±10 cm around the center of the interaction region. 

Primary charged particles are defined as all prompt particles produced in the collision, including de- 
cay products, except those from weak decays of strange hadrons. Selections based on the number of 
space points and the quality of the track fit, as well as on the distance of closest approach to the recon- 
structed vertex, are applied to the reconstructed tracks (see |8| for details). The efficiency and purity of 
the primary charged particle selection are estimated from a Monte Carlo simulation using the DPMJET 
event generator |[T2i with particle transport through the detector using GEANT3 |[T3l . The systematic 
uncertainties on corrections are estimated via a comparison to a Monte Carlo simulation using the HI- 
JING event generator lfT4ll . The overall primary charged particle reconstruction efficiency (the product 
of tracking efficiency and acceptance) for |T7iab| < 0.8 is 79% at pj = 0.5 GeV/c, reaches 81% at 0.8 
GeV/c and decreases to 72% for pj > 2 GeV/c. From Monte Carlo simulations it is estimated that the 
residual contamination from secondary particles is 1.6% at pj = 0.5 GeV/c and decreases to about 0.6% 
for pt>2 GeV/c. 

The transverse momentum of charged particles is determined from the track curvature in the magnetic 
field of 0.5 T. The pj resolution is estimated from the space-point residuals to the track fit and verified 
by the width of the invariant mass of Kg mesons reconstructed in their decay to two charged pions. For 
the selected tracks the relative pj resolution is 1.3% at pj = 0.5 GeV/c, has a minimum of 1.0% at 
Pt = I GeV/c, and increases linearly to 2.2% at pj = 20 GeV/c. The uncertainty on the pj resolution is 
±0.7% at Pj = 20 GeV/c, leading to a systematic uncertainty on the differential yield of up to 3% at this 
Pt value. 

Table 1: Systematic uncertainties on the /jj-differential yields in p-Pb and pp collisions for |77cms| < 0.3. The 
quoted ranges span the pj dependence of the uncertainties. 



Uncertainty 


Value 


Event selection 


1.0-2.0% 


Track selection 


0.9-2.7% 


Tracking efficiency 


3.0% 


Pj resolution 


0-3.0% 


Particle composition 


2.2-3.1% 


MC generator used for correction 


1.0% 


Secondary particle rejection 


0.4-1.1% 


Material budget 


0-0.5% 


Acceptance (conversion to T]cms) 


0-0.6% 


Total for p-Pb, pj-dependent 


5.2-5.5% 


Normalization p-Pb 


3.1% 


Total for pp, /^T-dependent 


7.7-8.2% 


Normalization pp 


3.6% 


Nuclear overlap (rpPb) 


3.6% 



Due to the different energy per nucleon of the two colliding beams, imposed by the two-in-one magnet 
design of the LHC, the nucleon-nucleon cms moves with a rapidity 3'nn = 0.465 in the direction of the 
proton beam. As a consequence, the detector coverage, |T7iab| < 0.8, implies, for the nucleon-nucleon 
cms, roughly —0.3 < T]cms < 1-3. The calculation of T]cn,s = Viab ±3'nn is accurate only for massless 
particles or at high pj. Consequently, the differential yield at low pj suffers from a distortion, which 
is estimated and corrected for based on the particle composition in the HIJING event generator. For 
Pj = 0.5 GeV/c, the correction is 1% for |T]cms| < 0.3 and reaches 3% for 0.8 < r/cms < 1-3. The 
systematic uncertainties were estimated by varying the relative particle abundances by factors of 2 around 
the nominal values. The uncertainty is sizable only at low pj and is dependent on Tjcms- It is 0.6% for 
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|T]cms| < 0.3, 4.3% for 0.3 < T]cms < 0.8, and 5.1% for 0.8 < T]cms < 1-3. 

The systematic uncertainties on the pj spectrum are summarized in Table [T] for Irjcmsl < 0.3. The total 
uncertainties exhibit a weak pj and T]cms dependence. The total systematic uncertainties range between 
5.2% and 5.5% for |T]cras| < 0.3 and reach between 5.6% and 7.1% for 0.8 < T]cms < 1-3. 

In order to quantify nuclear effects in p-Pb collisions, a comparison to a reference pj spectrum in pp 
collisions is needed. In the absence of a measurement at ^/s =5.02 TeV, the reference spectrum is ob- 
tained by interpolating or scaling data measured at =2.76 and 7 TeV. For pj < 5 GeV/c, the measured 
invariant cross section for charged particle production in pp collisions, d^a^^ /drjdpj, is interpolated bin- 
by-bin, assuming a power law dependence as a function of y^. For pj > 5 GeV/c, the measured data 
at ^/s = 7 TeV is scaled by a factor obtained from NLO pQCD calculations ifTSl . For pt < 5 GeV/c, 
the largest of the relative systematic uncertainties of the spectrum at 2.76 or 7 TeV is assigned as the 
systematic uncertainty at the interpolated energy. For pj > 5 GeV/c, the relative difference between 
the NLO-scaled spectrum for different choices of the renormalization /i« and factorization jXp scales 
(Hr = Hf = Pj, pj/2, 2pj) is added to the systematic uncertainties on the spectrum at 7 TeV. In ad- 
dition, an uncertainty of 2.2% is estimated comparing the interpolated and the NLO-scaled data. The 
total systematic uncertainty range from 7.7% to 8.2% for 0.5 < pj <20 GeV/c. The NLO-based scaling 
of the data at ^ = 2.76 TeV gives a result well within these uncertainties. More details can be found 
inHa. 




1 10 
Pj (GeV/c) 



Fig. 1: Transverse momentum distributions of charged particles in NSD p-Pb collisions for different pseudorapid- 
ity ranges (upper panel). The spectra are scaled by the factors indicated. The histogram represents the reference 
spectrum in pp collisions (see text). The lower panel shows the ratio of the spectra at forward pseudorapidities to 
that at |T7ems| < 0.3. The vertical bars (boxes) represent the statistical (systematic) errors. 

The final pp reference spectrum is the product of the interpolated invariant cross section and the average 
nuclear overlap (rpPb), calculated employing the Glauber model 11171 . which gives (rpPb) = 0.0983 ± 
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0.0035 mb ^ The uncertainty is obtained by varying the parameters in the Glauber model calculation, 
see ifTTl . 

The Pi spectra of charged particles measured in NSD p-Pb collisions at ^snn = 5.02 TeV are shown 
in Fig. [T] together with the interpolated pp reference spectrum. At high pj, the pj distributions in p-Pb 
collisions are similar to that in pp collisions, as expected in the absence of nuclear effects. There is an 
indication of a softening of the pj spectrum when going from central to forward pseudorapidity. This 
is a small effect, as seen in the ratios of the spectra for forward pseudorapidities to that at Irjcmsl < 0.3, 
shown in Fig. [T] (lower panel). Calculations with the DPMJET event generator |[T2l . which predict well 
the measured dA'ch/dT7iab IfTTl . overpredict the spectra by up to 33% for pj < 0.7 GeV/c and underpredict 
them by up to 50% for pj > 0.7 GeV/c. 

In order to quantify nuclear effects in p-Pb collisions, the pj -differential yield relative to the pp reference, 
the nuclear modification factor, is calculated as: 



dXfVdT7dpT 
(rpPb)d2c7PP/dr7d/.T' 



(1) 



where N^^^ is the charged particle yield in p-Pb collisions. The nuclear modification factor is unity for 
hard processes which are expected to exhibit binary collision scaling. For the region of several tens of 
GeV, binary collision scaling was experimentally confirmed in Pb-Pb collisions at the LHC by the recent 
measurements of observables which are not affected by hot QCD matter, direct photon [ 18], Z° [ 19], and 
[201 production. The present measurement in p-Pb collisions extends this important experimental 
verification down to the GeV scale and to hadronic observables. 
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Fig. 2: The nuclear modification factor of charged particles as a function of transverse momentum in NSD p-Pb 
collisions at y^snn = 5.02 TeV. The data for |Tjcms| < 0.3 are compared to measurements |8| in central (0-5% 
centrality) and peripheral (70-80%) Pb-Pb collisions at -^/snn = 2.76 TeV. The statistical errors are represented 
by vertical bars, the systematic errors by (filled) boxes around data points. The relative systematic uncertainties on 
the normalization are shown as boxes around unity near pj = for p-Pb (left box), peripheral Pb-Pb (middle box) 
and central Pb-Pb (right box). 

The measurement of the nuclear modification factor /?ppb for charged particles at | T]cms I < 0.3, is shown in 
Fig.|2l The uncertainties of the p-Pb and pp spectra are added in quadrature, separately for the statistical 
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Fig. 3: Transverse momentum dependence of the nuclear modification factor 7?ppb of charged particles measured in 
p-Pb collisions at y^s^ = 5 .02 TeV. The ALICE data in | Tjcms | < 0.3 (symbols) are compared to model calculations 
(bands or lines, see text for details; for HIJING, DHC stands for decoherent hard collisions). The vertical bars 
(boxes) show the statistical (systematic) errors. The relative systematic uncertainty on the normalization is shown 
as a box around unity near pj — 0. 



and systematic uncertainties. The total systematic uncertainty on the normalization, quadratic sum of the 
uncertainty on (rpPb), the normalization of the pp data and the normalization of the p-Pb data, amounts 
to 6.0%. 

In Fig. [2] we compare the measurement of the nuclear modification factor in p-Pb to that in central (0- 
5% centrality) and peripheral (70-80% centrality) Pb-Pb collisions at y^^NN = 2.76 TeV [8]. Rpp\, is 
consistent with unity for /?t ^ 2 GeV/c, demonstrating that the strong suppression observed in central 
Pb-Pb collisions at the LHC ||6]-[8l is not due to an initial-state effect, but rather a fingerprint of the hot 
matter created in collisions of heavy ions. 

The so-called Cronin effect 1211 (see ll22l for a review), namely a nuclear modification factor above unity 
at intermediate pj, was observed at lower energies in proton-nucleus collisions. In d-Au collisions at 
y^^NN = 200 GeV, /?dAu reached values of about 1.4 for charged hadrons in the pj range 3 to 5 GeV/c 
|[23l - l26l . The present measurement clearly indicates a smaller magnitude of the Cronin effect at the LHC; 
the data are even consistent with no enhancement within systematic uncertainties. 

Data in p-Pb are important also to provide constraints to models. For illustration, in Fig. [3] the mea- 
surement of RpPb at |T]cms| < 0.3 is compared to theoretical predictions. Note that the measurement is 
performed for NSD collisions. With the HIJING HH and DPMJET 1121 event generators, it is estimated 
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that the inclusion of single-diffractive events would lead to a decrease of /?ppb by 3-4%. Several pre- 
dictions based on the Colour Glass Condensate (CGC) model are available ||27] - |291 . The calculations 
of Tribedy and Venugopalan ETIl are shown for two implementations (rcBK and IP- Sat, see ETl for 
details). The calculations within IP-Sat are consistent with the data, while those within rcBK slightly 
underpredict the measurement. The prediction of Albacete et al. [28], for the rcBK Monte Carlo model, 
is consistent with the measurement within the rather large uncertainties of the model. The CGC calcu- 
lations of Rezaeian 1291 . not included in Fig. |3] are consistent with those of |[27ll28l . The shadowing 
calculations of Helenius et al. ll30l . performed at NLO with the EPS09s Parton Distribution Functions 
and DSS fragmentation functions describe the data well (the calculations ai^e for tt"). The predictions by 
Kang et al. |31 1, performed within a framework combining leading order pQCD and cold nuclear matter 
effects, show RpPb values below unity for Pt>6 GeV/c, which is not supported by the data. The predic- 
tion from the HIJING 2.1 model Il32l describes, with shadowing, the trend seen in the data, although it 
seems that, with the present shadowing parameter Sg, the model underpredicts the data. The comparisons 
in Fig. |3]clearly illustrate that the data are crucial for the theoretical understanding of cold nuclear matter 
as probed in p-Pb collisions at the LHC. 

In summary, we have reported measurements of the charged-particle pj spectra and nuclear modification 
factor in p-Pb collisions at ^/s^ = 5.02 TeV. The data, covering 0.5 < pj <20 GeV/c, show a nuclear 
modification factor consistent with unity for Pt>2 GeV/c. This measurement indicates that the strong 
suppression of hadron production at high pj observed at the LHC in Pb-Pb collisions is not due to an 
initial-state effect, but is the fingerprint of jet quenching in hot QCD matter. 
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